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Abstract 

The DIRAC spectrometer has been commissioned at CERN with the aim of de- 
tecting iT~^Tr~ atoms produced by a 24 GeV/c high intensity proton beam in thin 
foil targets. A challenging apparatus is required to cope with the high interaction 
rates involved, the triggering of pion pairs with very low relative momentum, and 
the measurement of the latter with resolution around 0.6 MeV/c. The general char- 
acteristics of the apparatus are explained and each part is described in some detail. 
The main features of the trigger system, data-acquisition, monitoring and setup 
performances are also given. 

Key words: DIRAC experiment, double arm spectrometer, pion scattering, 
experimental techniques, elementary atom 
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1 Introduction 



The DIRAC experiment aims to measure the ground state hfetime of 7r"'"7r~ 
atoms with 10% precision, using the 24 GeV/c proton beam of the CERN 
Proton Synchrotron. The atom Ufctimc is a consequence of the strong inter- 
action at low energy and it is determined by the charge exchange amphtude 
TT+TT" 7r°7r° very close to threshold. The probability of this process is pro- 
portional to the square of the difference of S-wave tttt scattering lengths with 

I 1 2 

isotopic spin and 2, |ao — 02! ■ The relation between the lifetime and |ao — 02! 
is model-independent [1]. The pion scattering lengths have been calculated in 
the framework of chiral perturbation theory with a precision of a few percent: 
ao = 0.220 ± 0.005 and aa = -0.0444 ± 0.0010. Using these values one can 
predict the pionium lifetime [2]: (2.9 ± 0.1) x 10~^^s. In order to determine 
|ao — 02! down to 5%, the lifetime has to be measured within 10% accuracy. 
Such a measurement would provide a crucial test for the understanding of 
chiral symmetry breaking in QCD. 

Pionium atoms (^27r) sue produced in proton-nucleus interactions. After pro- 
duction these relativistic atoms may either decay into 7r°7r° or get excited to 
higher quantum numbers, or break up (be ionised) in the target material where 
they are produced. In the case of break-up, characteristic pion pairs ("atomic" 
pairs) emerge. These pairs have a low relative momentum in their centre of 
mass system (Q < 3 MeV/c), very small opening angle (^ < 3 mrad) and 
nearly identical energies in the laboratory system. A high resolution magnetic 
spectrometer is then required [3] to split up the pairs and measure their rel- 
ative momentum with sufficient precision (0.6 MeV/c) to detect the pionium 
signal superimposed on the substantial background of "free" tt+tt" pairs pro- 
duced in inclusive proton-nucleus interactions. A previous experiment, using 
internal proton beam, has reported observation of pionium atoms [4]. 

The total number of produced 7r"'"7r~ atoms is related by an exact expression 
to the number of free pion pairs with low relative momenta. For a given target 
material and thickness the ratio of observed atomic pairs to the total number 
of produced atoms, i.e. the atom breakup probability, depends on the lifetime 
in a unique way [5]. 



2 General layout of the experimental setup 

The DIRAC experimental setup [6,7] is located at the T8 proton beam line 
of 24 GeV/c momentum in the East Hall of the PS accelerator at CERN. 
The isometric view of the setup is shown in Fig.l. The DIRAC apparatus 
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is designed to detect charged pion pairs with high resolution over the pair 
relative momentum. It became operational at the end of 1998 and has been 
collecting data since the middle of 1999. 

Fig. 1. Isometric view of the DIRAC setup. The radiation shielding boundaries are 
shown on the floor (each division marked on the boundary corresponds to 1 meter) . 

Fig. 2. Side view of the DIRAC setup. The secondary particle channel is inclined 
by 5.7° with respect to the primary proton beam. 

The setup consists of the proton beam line, target station, secondary particle 
vacuum channel, spectrometer magnet and detectors placed upstream and 
downstream the analysing magnet. Free and atomic 7r"'"7r~ pairs produced in 
the target enter the secondary particle channel which is tilted upwards by 5.7° 
with respect to the proton beam (Fig. 2). At the end of the secondary particle 
channel the spectrometer magnet is installed, also tilted by 5.7° together with 
all the downstream detectors. 

The top view of the setup is shown in Fig. 3. The upstream section of the 
secondary particle channel between the target station and the spectrometer 
magnet is instrumented with the following detectors: microstrip gas chambers 
(GEM/MSGC), scintillating fibre detector (SFD) and scintillation ionisation 
hodoscope (IH). 

Downstream the spectrometer magnet the setup splits into two identical arms 
for detection and identification of positive and negative charged particles. The 
angle between each arm and the spectrometer symmetry axis is 19°. Along 
each arm the following detectors are located: drift chamber system (DC), 
vertical scintillation hodoscope (VH), horizontal scintillation hodoscope (HH), 
gas Cherenkov counter (CH), preshower detector (PSH) and muon detector 
(MU). 



3 Beam lines, spectrometer magnet and radi- 
ation shielding 

3.1 Proton beam and target station 



To extract protons from the PS to the T8 beam line a slow ejection mode is 
used. The beam is extracted in spills of ~ 400-500 ms duration. During data 
taking, between 1 to 5 cycles per PS super-cycle of 14.4 19.2 s duration are 
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Fig. 3. Schematic top view of the DIRAC spectrometer. Moving from the tar- 
get station towards tlic magnet there arc 4 planes of microstrip gas chambers 
(GEM/MSGC), 3 planes of scintillating fibre detectors (SFD) and 4 planes of ionisa- 
tion hodoscope (IH). Downstream the dipole magnet, on each arm of the spectrome- 
ter, there are 4 stations of drift chambers (DC), vertical and horizontal scintillation 
hodoscopes (VH, HH), gas Cherenkov counter (CH), preshower detector (PSH) and, 
behind the iron absorber, muon detector (MU). 

delivered to DIRAC. The proton beam intensity was set to (0.6 — 1.0) • 10^^ 
protons per spill, depending on the target used. 

The PS proton beam line includes two bending magnets deflecting the beam at 
an angle of 76 mrad towards the final straight T8 section, corrector magnets 
performing horizontal and vertical steering and quadrupole magnets which 
focus the beam on the experiment target. The dimension of the beam spot 
at the target location are x = 1.6 mm, y = 3.2 mm at 2a level. The diver- 
gence of the beam is about 1 mrad. The nominal momentum of the extracted 
beam is 24 GeV/c, with instantaneous momentum spread close to 0.08% at 
2a. The design of the proton beam line optics has been optimised using the 
TRANSPORT simulation code [8] by the CERN PS Division. 

Downstream the target the proton beam travels in a vacuum channel below 
the spectrometer magnet and detectors and finally is absorbed by a beam 
dump. 

To measure the beam intensity and to tune the beam position on the target, 
the beam line is equipped at several locations with secondary emission cham- 
bers and luminescent screens with TV cameras. One additional beam position 
detector (centroid) [9] is installed close to the target station. 



The target station houses a remote controlled device with 12 holders for the 
targets, including an empty holder and a luminescence screen. During data 
taking, targets made of Pt (28 /xm thick), Ni (94 /xm and 98 /xm thick) and 
Ti (250 yum thick) were used. 

The DIRAC experiment is sensitive to particles outside the beam core (halo) , 
because the target is very thin (nuclear target efficiency is < 10~^), and the 
upstream detectors are placed very close (18-26 cm in the vertical direction) 
to the primary proton beam. The halo is originated from scattering of primary 
protons on the sphtter blades, and a special optics has been designed to de- 
crease the background halo to a neghgible level. The ratio of detector counting 
rates with the target in place to those with an empty holder was measured to 
be ~25. 



3.2 Secondary particle channel and spectrom- 
eter magnet 

The secondary particle channel [6] [7] is placed at an angle of 5.7° relative 
to the proton beam and consists of two vacuum volumes, as shown in Fig. 3. 
The first one is a 2 m long, 611 mm diameter tube, located immediately 
downstream the target station, common to both the proton beam line and 
the secondary particle channel. Secondary particles exit this tube through 
a 200 mm diameter window, made of 250 /xm thick mylar film. The second 
volume, located at ~3.5 m from the target, consists of a cylindrical vacuum 
section, containing a collimator, attached to a 2.7 m long fiat vacuum chamber 
placed between the spectrometer magnet poles. A 1.5 m long air gap between 
the two vacuum sections allows insertion of the upstream detectors. The angu- 
lar aperture of the secondary particle channel is determined by the collimator 
and is equal to ±1° in horizontal and vertical directions resulting in a solid 
angle acceptance of 1.2 ■ 10^'^ sr. The fiat chamber is ended with a 0.68 mm 
thick Al outlet window of 2.0 x 0.4 m^ dimensions {W x H). 

The spectrometer dipole magnet (magnetic field B=1.65 T, field integral BL = 
2.2 T-m) has an aperture of 1.55 x 0.50 m^ {W x if). To reduce the stray field, 
two magnetic screens are fixed near its entrance and exit. 

Table 1 summarises the material thicknesses (in units of radiation length 
xlO""^) encountered by secondary particles before they reach the DC system 
where their momenta are measured. 
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Table 1 

Material contributions along the secondary particle channel in units of radiation 
length xlO~^. 



Ni-target 


33.5 


Mylar window 


8.7 


4 olanes GEM/MSGC 


224.1 


3 planes SFD 


260.0 


4 planes IH 


153.1 


air gap 


34.7 


Mylar window 


8.7 


Al-window 


76.4 


Total 


765.7 



3.3 Beam dump and radiation shielding 

Neutron and gamma fluxes may cause a serious problem to sensitive elements 
of the setup. To estimate their effect a simulation of the background radiation 
flux in the full experimental apparatus has been performed [10]. The results 
were used to optimise the design of the radiation shielding which is shown in 
Fig. 4. 




I DIRAC I 



Fig. 4. The DIRAC setup on the T8 PS extraction line and the radiation shielding. 

At the end of the T8 beam line primary protons are absorbed by an iron beam 
dump. To decrease the background gamma and neutron fluxes from the beam 
dump towards the detectors, a dedicated radiation shielding has been adopted. 
It includes installation of a graphite core into the beam dump concrete 
wall near the beam dump and, at 3 m distance, another iron-concrete-iron 
wall, both with holes to allow passage of the proton beam pipe. 

The downstream detectors are shielded in addition from background secondary 
particles produced on the primary proton pipe and surrounding elements. For 
this purpose aim thick iron wall is installed between the upstream detector 
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region and the spectrometer magnet. In addition, collimators are inserted both 
in the primary proton beam pipe and in the secondary particle channel. The 
presence of the collimator in the proton beam line determines a reduction of 
the background rate by a factor of 2. 

A radiation shielding encloses the whole DIRAC experimental apparatus to 
protect the surrounding East Hall area from irradiation. It has been designed 
according to the maximum flux of 2.7 x 10^° incident protons per second, 
in respect of the CERN safety regulations. Being the apparatus located in a 
fully enclosed area, it has become necessary to provide the experimental area 
with cooling and ventilation equipment to prevent overheating of detectors 
and electronics. 



4 Large Detectors upstream the magnet 



Two tracking devices have been installed in the secondary particle channel 
before the magnet and the collimator : the GEM/MSGC and the SFD. They 
are used to improve the resolution on the measurement of the longitudinal and 
transverse components of the relative momentum of pion pairs as determined 
by the drift chambers tracking system and the nominal position of the beam 
at the target center. 

At the same time these detectors allow to select particle pairs originated by 
primary interactions at the target from the background of secondary inter- 
actions and particle decays. The ensemble of MSGC/GEM+SFD constitutes 
a tracking system with 7 detector planes with 2 stereo angles which provides 
adequate space resolution to reach the limit of multiple scattering in the target 
material. 

To increase the detection capability on close-lying tracks, an lonisation Ho- 
doscope (IH) is installed downstream the SFD, with the purpose of detecting 
pion pairs with a too small opening angle to be resolved by the tracking detec- 
tors. This is achieved by a detailed pulse-hight analysis of the double ionisation 
produced by the particle pairs in 4 layers of scintillation counters. 

A general picture of the above-mentioned detectors, as they are installed be- 
tween the first vacuum chamber and the secondary particle channel, can be 
seen in Figure 5. 
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Fig. 5. Photography of the three detectors installed upstream the magnet, between 
the first vacuum chamber (right-hand side) and the secondary particle channel. 
From right to left, the GEM/MSGC, SFD and IH detectors can be found. The 
primary proton beam line can be appreciated at the bottom. 

5 The GEM/MSGC detector 

This detector performs particle tracking at a distance of 2.4 m from the inter- 
action point. It is a proportional gas detector, based on the principle of the 

Gas Electron Amplifier (GEM) [11] [12], complemented with a second ampli- 
fication and readout stage provided by Micro Strip Gas Chambers (MSGC) 
[13] [14]. A more complete description of the detector and its performance is 
being prepared in a separate publication [15]. 

It measures particle coordinates in 4 planes along the direction of the incoming 
particle: X,Y,U,V, with orientations 0, 90, 5, 85 degrees, respectively, where 
the degrees are defined by microstrips running vertically (X-coordinate) . The 
stereo angles allow resolution of ghost combinations for two or more particles. 
With a single-hit space resolution close to 54 /im, this detector provides a 
precise measurement of the pion pair angular opening, ultimately limited by 
multiple scattering in the thin target. 



5.1 Detector concept 

Each chamber has active area 10.24x10.24 cm^, and consists of a drift elec- 
trode, a GEM foil and a MSGC sensor. The GEM plane is evenly spaced from 
the other two with a uniform gap of 3 mm, as indicated in Fig. 6. The drift 
electrode is made of a Chromium-coated thin glass (200 /xm). The GEM is a 
50 iim thick kapton foil copper-cladded on both sides with a 4 /xm thick Cu 
layer. The etching pattern is characterised by 50 /im wide holes, 140 /xm apart 
[16]. Application of a potential difference of 400 V between the two metal 
layers {Vi = —1600 V, V2 — —2000 V) produces electron amplification by a 
factor of 30 [12]. 

The MSGC sensor consists of 200 //m pitch alternating Chromium strips, with 
9 /im and 100 /xm anode and cathode width, respectively. They are implanted 

on a bare DESAG D263 substrate. Applied voltages are: —410 V on cathodes, 
and —3000 V on the drift electrode, whereas anode strips are set to ground. 
The gas employed is a mixture of Ar-DME (60/40). Under these conservative 
conditions an overall detector gain of approximately 3000 is achieved. 
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DESAG D263 200 \m glass with 
.2 nm Chromium 




Cr microstrips: 9[im anodes and 100 |im cathodes Jim copper) 

200 um Ditch: substrate D263 300 um 



Fig. 6. Layout of the GEM/MSGC detector. 

5.2 Readout electronics 



It is composed of three parts: the analog boards, the control boards and the 
VME acquisition modules [17]. 

The analog boards contain the front-end electronics, and also serve as a me- 
chanical support for the detector. They host 16 Analog Pipeline Chips (APC) 
[18], bonded on a multi-layer hybrid to the detector fan-out strips. Each APC 
has 64 inputs, with a charge sensitive preamplifier followed by analog pipeline 
of 32 capacitors, which is run at 10 MHz. In order to accelerate the serial 
readout, only 32 inputs are connected to the same chip. The analog pulse is 
formed by charge subtraction from consecutive capacitors, at positions given 
by the delayed trigger signal. This robust procedure, which allows easy inte- 
gration of a large number of channels, limits however the time resolution of 
the readout to a window of approximately 200 ns. 

The control boards store the digital sequence needed for APC control and 
perform fast digital conversion (AD9048) of the 16 series of analog signals. 
These two functions are controlled by FPGA (Fast Programmable Gate Array) 
incorporated onto the board. 

The VME modules handle the trigger signal and perform real time pedestal 
subtraction and zero suppression. One module is divided into 4 segments, each 
controlled by an independent FPGA. A Fast Clear function is implemented to 
be used in connection with the highest level of trigger (see trigger chapter), 
that stops the readout of current event and enables new triggers. 
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Fig. 7. a) GEM/MSGC strip multiplicity per hit-cluster b) cluster multiplicity per 
event in one plane (X). The dotted line shows all clusters registered by the data 
acquisition, whereas the continuous line shows only those having a time tag, es- 
tablished when the cluster has a corresponding hit in the SFD aligned with the 
interaction point. 

5.3 Performance 



Neighbouring hit strips are pattern-recognised as clusters. The hit multiplicity 
in them (a) and the cluster multiplicity per plane (b) are shown in Fig. 7. 
Single hit resolution was determined during the commissioning run in April 
2000 by setting all planes parallel to each other, and 54 yum is a typical value, 
as shown in Fig. 8. Two detector planes were installed in 1999 and the full 
set in 2000. The number of detector dead channels is around 1% and no 
significant deterioration due to radiation has been observed since then. The 
average efficiency for a standard detector is 93%. When at least 4 signals 
among 7 detectors are required to make a track, this is sufficient to provide 
99% overall tracking efficiency upstream the magnet. 



6 The Scintillating Fibre Detector 



The Scintillating Fibre Detector (SFD), together with the GEM/MSGC detec- 
tor, enables particle tracking to be performed upstream the magnet, with the 
required space and time resolution. In addition, it provides topological trig- 
ger capabilities [19] for rejection of pairs with relative distances larger than 
9 mm at the detector location (this feature was used at the early stage of the 
experiment, see a footnote in section "Trigger system"). 
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Fig. 8. Space resolution of the GEM/MSGC detector measured in a dedicated 
beam-test with 4 planes parallel to each other. 



6.1 Detector concept 



The SFD consists of three fibre planes to measure the X-, Y- and U-coordinates 
of incident particles ^ . The SFD covers a 105x105 mm^ area, each scintillating 
fibre (SciFi) array consists of several layers of KURARAY fibres. Five (for 
planes X and Y) or three (for plane U) fibres forming one sensitive column 
are mapped onto one channel of position-sensitive photomultiplier (PSPM). 
A layout illustrating the main characteristics of the SFD (for a plane with 5 
fibres per column) and PSPM mapping is shown in Fig. 9. 



^ The SFD operational in DIRAC until the end of 2001 run consisted of X and Y 
planes only. The U-plane described here was installed in early 2002. 
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Fiber array Fiber position on PSPIVi 



Fig. 9. The SFD principal structure. A 16-channel fragment is shown. 

The fibre columns pitch is 0.44 mm, thus allowing some overlap to optimise 
efficiency. Each SciFi is connected by optical epoxy to a clear fibre fight guide 
of ~300 mm length. The edges of clear fibres, assembled in a bundle, are 
carefully polished and glued into the holes of a square black plate which fixes 
the fibre positions on the PSPM photocathode; no optical grease is used. The 
far end of the SciFi array is tightly connected to a mirror made of aluminised 
mylar. With this arrangement the fight attenuation along the SciFi was found 
to be negfigible. The details of the fibre arrays are fisted in Table 2. 

Table 2 

SFD planes specification. 





Fibre 


Fibre 


Length 


Fibres/ 


Xo 


Num. of 


Num. of 




type 


$ (mm) 


mm 


column 


% 


channels 


PSPM 


X 


SCSF38 


0.50 


130 


5 


0.8 


240 


15 


Y 


SCSF38 


0.50 


130 


5 


0.8 


240 


15 


U 


SCSF78M 


0.57 


150, 


3 


1.0 


320 


20 








130, 70 











The SFD individual planes are installed close to each other. Direction of fibres 
is orthogonal in planes X and Y while plane U is rotated by 45 degrees. The 
U-plane has 5 sections of different fibre lengths to cover roughly the same area 
like the X- and Y-planes. 

A 16-channel metal dynode position-sensitive photomultiplier tube, Hama- 
matsu H6568, has been selected as photosensor. This photomultiplier is char- 
acterised by good timing properties (rise time ~0.7 ns), low noise (1-2 pulses/s 
at a nominal detection threshold) and perfect single photoelectron spectrum. 
It has been modified to monitor the amplitude from the last dynode for cal- 
ibration purposes. The level of optical cross talk among the PSPM channels 
was found to be ~1% (with a 1.6 mm diameter light spot on the photocath- 
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ode). The measured linear range of this tube (for the hnear bleeder) extends 
up to 15 photoelectrons at 950 V. 

The detector concept with almost independent channels was developed in the 
context of the RD-17 project [20]. The results from the SFD study with a test 
beam of low intensity, reported in [21], are the following: 



— light output 6-10 photoelectrons, 

— average detection efficiency 98.4%, 

— r.m.s. of the detection efficiency 1%, 

— average hit multiplicity 1.1, 

— spatial resolution (cr) 127 /xm, 

— time resolution (a) 0.65 ns. 




SciFi TDC spectrum for Y140 


Chi2/ndf = 121.1 /92 
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Fig. 10. SFD raw time spectra for e'^e" trigger data. Here and in Fig. 11 the 
horizontal scale is in TDC channels, the channel width is 0.5 ns. 
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6.2 Readout electronics 



A dedicated electronic circuit (PSC) has been custom developed to provide 
signal discrimination with dynamic rejection of cross-talk in adjacent channels 
using the peak-sensing technique. A detailed description of a 32-channel PSC 
module is given in [22]. Discrimination of a channel is given by the condition 
2Ai — Ai_i — Aj+i > Athr, where Ai are channel signal amplitudes and Af^r 
defines the threshold value. The PSC algorithm provides efficient detection 
of double tracks from time correlated particle pairs (up to ~5 ns time dif- 
ference) when the relative distance between the two tracks is larger than the 
fibre column pitch. However, when adjacent fibre columns are crossed by two 
particles simultaneously, then the PSC algorithm leads to a suppression (by 
20-40%) of the detected yield of double track events. In such cases a signal is 
detected for only one of two hit columns. The detection efficiency depends on 
the detector light output in a specific channel. For particle pairs with relative 
time differences greater than 10 ns (accidental pairs) the PSC behaves as an 
ordinary leading edge discriminator. 



SciFi TDC spectrum for X131 




Chi2/ndf= 62.29/37 
Level =40.83+- 1.216 
Constant = 166.2 +- 8.895 
Mean = 60.52 +- 0.1038 
Sigma = 2.251 +- 0.1193 



90 
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TDC spectrum for Y1 40^ chi2 / ndf = 62.45 / 32 

Level =42.76+- 1.344 
Constant = 126 +- 8.073 
Mean = 54.54 +- 0.1284 
Sigma = 2.124 +- 0.1367 




Fig. 11. SFD raw time spectra for tt+tt trigger data. 

The front-end electronics is located inside the beam area close to the detector. 
The PSC output ECL signals are sent to LeCroy 3377 multi-hit TDC. This 
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solution does not require additional delay lines to adjust the timing of the 
SFD with respect to the trigger timing. 



6.3 Performance in the experiment 

Due to the high flux of particles at the position of the detector close to the 
target, and to the presence of a non-negligible amount of inclined tracks associ- 
ated with secondary interactions in the channel, the SFD performances in the 
experiment slightly differ from the above mentioned. The detection efficiency 
is still high (around 98%), but the average hit multiplicity is near 5 in the 
50 ns time window of TDC (at a nominal beam intensity of 10^^ protons per 
spill impinging on a 94 ^m Ni target). The raw time spectra, obtained from 
e+e" and 7t~^7t~ events, are shown in Fig. 10 and Fig. 11, respectively, for two 
arbitrary SFD channels. The width of the distributions is dominated by the 
time jitter of the trigger signal. After off-line deconvolution of the trigger time 
jitter the resolution of the SFD is found to be a—O.S ns. 



7 The lonisation Hodoscope 

Charged pions originated from pionium breakup cross the upstream detectors 
at rather small relative distances. When the distance is less than the double 
track resolution of the upstream tracking devices, then only one hit is de- 
tected, thus making the event reconstruction ambiguous. That is why another 
technique based on a measurement of the ionisation loss is used as well. 

A dedicated lonisation Hodoscope (IH) [23] has been built to separate double 
ionisation signals produced by close pion pairs incident on the same scintillat- 
ing slab, from single ionisation signals produced by one particle. In this way, 
the uncertainties resulting from the inefficiency in detecting two tracks with 
relative distance approaching zero can be significantly reduced. 

The lonisation Hodoscope described here was installed in 2001 to replace a 
previous version of a similar detector type [24], consisting of only two planes 
with 16 slabs of 2 mm thickness, oriented in the vertical direction. 

The present IH detector is a scintillation hodoscope consisting of 4 planes of 

11x11 cm^ sensitive area placed normally to the axis of the setup (Fig. 12), 3 m 
downstream the target. Two planes have vertically oriented slabs (planes X-A 
and X-B) whereas the other two have horizontal slabs (planes Y-A and Y-B). 
They are arranged in the following sequence, moving along the beam direction: 
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Fig. 12. Design of the IH scintillation plane. 1 - scintillators, 2 - light-guides, 3 - PM 
photocathodes. 




Fig. 13. Isometric view of the lonisation Hodoscope. 1 - scintillators, 2 - light-guides, 
3 - photomultipliers with shielding. 

X-A, Y-A, X-B, Y-B. This ordering has been chosen to minimise possible cross- 
correlations between signals in the planes (e.g. due to 5-electrons). Each plane 
is assembled from 16 plastic scintillating slabs made of fast scintillator (BC- 
408). Planes with the same slab orientation are shifted by a half-slab- width 
with respect to each other. The slabs are 11 cm long, 7 mm wide and 1 mm 
thick. They are connected to the PM photocathodes via 2 mm thick and 7 mm 
wide Incite hght guides (fig. 13). 

The front and rear surfaces of a slab are covered by a millipore film [25] for 
efficient hght collection. At the lateral surface of the slab, hght is reflected by 
a thin (30 ^m) aluminised black mylar film, which is used instead of millipore 
film in order to minimise the gaps between adjacent slabs. A typical gap 
between two adjacent slabs in this configuration is less than 70 fim. wide. 

Scintillation light is detected by FEU-85 photomultiphers with 25 mm diame- 
ter photocathodes. Photomultipliers are assembled by 16 units into a compact 
set, allowing independent replacement of each PM. Photocathodes are in op- 
tical contact with the wide side of a light guide instead of the traditional 
butt-end readout. This improves the hght collection efficiency by about 50%. 



17 




400 



ADC channels 



Fig. 14. Typical ADC spectra for single (solid line) and double (dashed line) ioni- 
sation loss from particles crossing one IH scintillating slab. 




Fig. 15. Contamination of single ionisation amplitudes as a function of losses of 
double ionisation as obtained from the analysis of the spectra of Fig. 14. 



As the detector is highly loaded by intense particle flux, the last 4 PM-dynodes 
are fed by an additional power supply to ensure a constant PM amplification 
throughout the spill. 

Signal amplitude and time are digitised by LeCroy ADC 4300B and TDC 3377 
modules, respectively. The time resolution of the IH detector is better than 

1 ns. The typical response of one IH channel to close particle pairs incident on 
one scintillating slab and to single particle is shown in Fig. 14. If a threshold is 
set to retain 90% of the double ionisation signal from pairs, the contamination 
from single particle amplitudes is less than 15% (Fig. 15). 
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8 Drift Chambers 



8.1 General layout and characteristics 

The drift chamber system is used to perform particle tracking downstream 
the dipole magnet. The system is designed to sustain a high particle fluency 
in the secondary channel, reaching 10 kHz/cm^ at the innermost region. 

A two-arm solution has been chosen, except for the first chamber which is a 
single large module (DC-1) designed with two separated sensitive areas 0.8 x 
0.4 m^ each. This chamber provides 6 successive measurements of the particle 
trajectory along the coordinates X,Y,W,X,Y,W, where W is a stereo angle 
with inclination 11.3° with respect to the X-coordinate. DC-1 is instrumented 
with 800 electronic channels. 

Each of the two arms consists of 3 chamber modules, of identical design, mea- 
suring coordinates X,Y (DC-2), X,Y (DC-3) and X,Y,X,Y (DC-4) following 
the direction of the outgoing particle. Their dimensions are 0.8 x 0.4 m^ (DC- 
2), 1.12x0.4 m^ (DC-3), and 1.28x0.4 m^ (DC-4). Both arms together contain 
1216 electronic channels. 

The distance between the center of the first half of DCl and the center of DC4 
provides a lever- arm of 1.6 m along the average particle path, having uniform 
spacing of chambers DC-2 and DC-3 along this path. Characteristics of the 
drift chamber system are summarised in Table 3. 



8.2 Drift chamber electrodes 

A schematic drawing of the sensitive element is shown in Fig. 16. The anode 
wires pitch is 10 mm, the distance L between the anode and cathode planes is 
5 mm. The cathode planes and potential wires are at equal voltages. As seen 
in the figure, a sensitive area, corresponding to each anode wire and limited 
by the cathode planes and potential wires, has a square (10 x 10 mm^) shape. 
In this case, with a suitable gas mixture, it is possible to achieve a linear 
behaviour of the drift function, except in a small region near the potential 
wire. 

Cathode planes are made of 20 thick carbon-coated mylar foils with a 
surface resistivity of about 400 Q per square. Such cathode foils provide stable 
chamber operation due to a high work function of the carbon coating and, 
being thin, add only small amount of material along the particle path. 
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Table 3 

General properties of the DC modules. 



Module 


Sensitive 


Measured 


Number of 


type 


area, cm^ 


coordinate 


planes 


DC-1 


40 X 80 


X 


2 




left arm 


Y 


2 






W 


2 




40 X 80 


X 


2 




right arm 


Y 


2 






W 


2 


DC-2 


40 X 80 


X 


1 
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DC-3 


40 X 112 
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1 
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Fig. 16. Schematic view of the wire chamber electrodes: AW - anode wires, PW 
potential wires, C - cathode foils. Dimensions are in mm. 



Anode and potential wires of 50 jim. and 100 jim diameter, respectively, are 
made of a copper-beryllium alloy. The rather large diameter of the anode wires 
has been chosen in order to operate the chambers at high current avalanche 
amplification mode. 
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8.3 Chamber design 



The chamber design is shown in Fig. 17 for the case of the DC-2 module. The 
module is a stack of aluminium and fibreglass frames, each of 5 mm thick- 
ness, fixed by screws. The fibreglass frames are the supports for the chamber 
electrodes (anode and potential wires and cathode foils). The two outer alu- 
minium frames in the stack are used to fix the mylar window, and the inner 
ones are the spacers between the fibreglass frames. Rigidity of the module is 
enforced by aluminium rectangular tubes screwed to the surface of the frame 
package. 
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Fig. 17. Design of the DC-2 module. Upper figure: general view. Lower figure: 
structure of the frame stack; X - X-plane, Y - Y-plane, C - cathode foils. 

Gas tightness of the chamber module is provided by rubber o-rings glued 
along the inner edges of the aluminium frames. Within a module, gas flows 
sequentially in the sub-volumes defined between cathode foils, by means of 
holes drilled on opposite sides of the fibreglass frames. 

The design of module DC-1 differs from the one shown in Fig. 17. The main 
difference, illustrated in Fig. 18, consists in the fact that DC-1 comprises, in 
a single gas volume, two sets of sensitive planes, placed symmetrically to the 
left and right hand side of the spectrometer axis. The middle zone, which is 
strongly irradiated by particles (mostly fast protons from target fragmenta- 
tion) , is made insensitive to the particle flux. The limiting edge of the sensitive 
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zones, close to the axis, can be varied. This is possible by means of a stripped 
structure of the neighbouring cathodes, which allows stepwise application of 
voltage. This design of the DC-1 module ensures little amount of material, by 
avoiding frames in the small angle region. 

X-plane 

































i 







Y-plane 



V///M 



ff-plane 



Fig. 18. Schematic view of the DC-\ module. Dotted areas show the sensitive 
regions of the Y- and VF-planes. Hatched areas mark the zones of the cathode 
strips which allow to change the width of the insensitive area in the central region. 



8.4 Chamber operation and performance 

The drift chambers operate in a high current avalanche mode. This mode is 
characterised by high pulse amplitude (about 1 mA), small pulse width (20 ns), 
and stable operation due to an efficiency plateau larger than 1 kV. The single 
hit efficiency is above 96% when the particle flux is about 10 kHz/cm^. The 
employed gas mixture is Ar(~ 50%) +iC4i?io(~ 50%) + i72O(0.5%), and the 
chamber operation voltage is 3.85 kV. 

A space-to-time relationship was extracted from the time spectrum and its 
integral distribution shown in Fig. 19, for a sample of clean events with a small 
amount of background hits. The integral distribution has been parameterised 
by a second order polynomial of the type: 

Z = ai X a2 X i*^. 

In this formula t* — txDC ~ to ~ ^t, where St is the signal propagation time 
along the anode wire. 

Study of the drift function parameters for different chamber planes at different 
beam intensities shows good stability of the above relation. For this reason the 
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same drift function parameters {cii and 02) were used for all chamber planes 
during the off-line track reconstruction procedure. 
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Fig. 19. Distribution of the drift time (upper) and its integral spectrum (lower) for 
the X4-plane. Horizontal scale is in TDC channels, bin width is 0.5 ns. 

Coordinate resolution of the DC system is illustrated in Fig. 20, where the 
distribution of differences between the predicted position and measured co- 
ordinates in one of the planes is shown (X4- plane right arm). The measured 
standard deviation, cr = 100 fim, is defined not only by the intrinsic chamber 
plane resolution, but also by the accuracy of the predicted track coordinates. 
Taking the latter into account the measured intrinsic space resolution of one 
plane is better than 90 /im. 

Tracking efficiency of the drift chamber system as a whole is about 99%, due 
to the fact that the requested number of hits per reconstructed track is less 
than the total number of sensitive planes crossed by a particle. 



8.5 Readout electronics 

The readout electronics of the drift chambers, which is a custom-made system 
[26] , provides data readout into the data collection memories and input to the 
trigger processor (see section "Trigger system" ) . 

The sensitive wire signals are digitised in the 16-channel multi-hit time-to- 
digital converter boards (TDC), which are plugged in the connectors mounted 
onto the chamber frames. This solution results in reduced number of electronic 
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Fig. 20. Distribution of differences between the measured and predicted 
X-coordinate for one drift chamber plane (XA). 



units, small number of cables and high noise immunity. The detection thresh- 
old in the TDC board can vary from 0.05 to 2 mA, the maximum number of 
hits per channel is 16. Least count of TDC is 0.5 ns, which corresponds to a 
drift distance of 25 /xm at a constant drift velocity of 50 /im/ns, well below 
the chamber intrinsic space resolution. 

A complete readout chain consists of TDC boards, bus drivers, readout con- 
troller and VME memory. Up to 8 TDC boards can be connected to the bus 
driver via common data and control buses, forming a segment of the DC read- 
out system. Similarly, up to 8 segments are connected to a readout controller 
forming a readout branch. The data of an event are stored in local data buffers 
until the higher level trigger decision is issued. If the event is accepted, the 
data arc serially transferred via the readout controller to the VME memories. 
To read out all the DC data 3 readout branches are used. The accepted data 
are transferred to the VME buffers within 5 us on average. The DC readout 
is fast enough compared to the global readout time of the experiment. 

For trigger purposes the DC readout system is equipped with fast data ports 
which directly transmit the hit wire numbers to a trigger track processor. This 
provides a minimum access time to the data, thus reducing the latency of the 
trigger system. Transmission to the track processor is performed in parallel. 
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9 The Vertical Hodoscopes (Time-of-Flight De- 
tector) 

The vertical hodoscope (VH) consists of an array of vertical scintillating slabs 
placed downstream the DC system. The VH system, together with the hori- 
zontal hodoscope, provides fast coincidence signals between the spectrometer 
arms necessary for the first level trigger. It is also used, in correlation with 
other detectors, in the definition of dedicated triggers for calibration purposes 
and of a higher level trigger for the selection of low Q events (see trigger 
section later). 
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Fig. 21. Time difference between charged positive and negative particles obtained 
from standard hadron trigger data. The central peak has a gaussian width of 193 ps, 
and the shaded area represents a 2a cut used to select prompt tttt events. The fiat 
background is originated from accidental pairs, not belonging to the same beam 
interaction. Note the shoulder on the right-hand side of the peak, due to ■K~p prompt 
pairs. 

A key function of this detector, which motivated a special design, is to provide 
a very accurate time definition of pion pairs originated from the same proton 
interaction (prompt pairs), in order to perform a clean separation (in off-line 
analysis) with respect to pairs in which the pions are produced at different 
times (accidental pairs). 
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Used as a time-of-flight detector it allows to identify p'K~ pairs in prompt 
events, as they might constitute a significant source of background to the 
7r"'"7r~ signal. 

The VH detector consists of two identical telescopes matching the acceptance 
of the drift chamber system previously described. Each telescope contains an 
array of 18 vertical scintillation counters. The scintillating material is BI- 
CRON BC420 and the slab dimensions are 40 cm length, 7 cm width and 
2.2 cm thickness. Scintillation light is collected at both ends by two 12-dynode 
Hamamatsu R1828-01 photomultipliers (chosen because of their small tran- 
sition time spread) coupled to fish-tail light guides. Voltage dividers were 
designed to provide sufficient bleeder current at particle rates up to 2x10^ 
Hz, while preventing any degradation of the time resolution. 

The front-end electronics was designed to minimise the time jitter. This is 
achieved by using LeCroy L3420 constant fraction discriminators, followed by 

CAEN C561 meantimers to provide a position independent time measurement. 
Hodoscope signals are delayed and transferred to the data acquisition system 
using shielded twisted-pair cables. Time digitisation is performed by LeCroy 
4303 time-to-FERA Converter followed by ADC 4300B. The least count of 
this ensemble is 62 ps. 

The VH single-hit detection efficiency is 99.5% for the positive, and 98.8% for 
the negative hodoscope arms. In Fig. 21 the distribution of the time difference 
between positive and negative pions in the spectrometer is shown. The ob- 
served ratio between prompt and accidental pairs in the 2a cut region around 
the peak is about 16. 

The overall time resolution of the system has been measured with e'^e~ pairs 
to be 127 ps per counter [27], which corresponds to 174 ps accuracy for the 
time difference between positive and negative arms (time-of- flight resolution) . 
The latter is shown in Fig. 22. 

The dedicated e+e~ calibration trigger selects e'^e" pairs from 7 conversions 
and Dalitz decays of 7r° which are almost synchronous in time, as the time of 
flight of e+e" pairs is momentum independent in the available setup range of 
momenta. 

This timing capability allows to separate tt+tt^ from n^p pairs in the mo- 
mentum range from 1 to 5 GeV/c, and from ttK pairs in the range from 1 to 
2.5 GeV/c, as illustrated in Fig. 23. 
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Fig. 22. Time difference spectrum for e^e~ pairs detected by the vertical ho- 
doscopes, after path length correction. Data come from a sample of e"'"e~ calibration 
triggers. 




Fig. 23. Correlation between the measured momentum of the positive particle and 
the VH time difference between the positive and negative spectrometer arm, taking 
into account the correction for the difference in path length. The accumulation 
bands correspond to 7r^7r+ (vertical band) and 7r~p (curved band) pairs. A small 
cluster of tt~K~^ pairs is also visible in the intermediate region. 
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Fig. 24. General layout of the horizontal hodoscopes. 

10 The Horizontal Hodoscopes 



The horizontal hodoscopc system (HH) is also separated into two arms, each 
covering an area of 40 x 130 cm^. The HH participates, together with the ver- 
tical hodoscope, to the definition of the first level trigger. Its response is used 
to apply a coplanarity criterion to track pairs hitting both detector arms. This 
trigger requirement selects oppositely charged particles with relative vertical 
displacement. At/, less than 7.5 cm. 

Each hodoscope consists of 16 horizontal extruded scintillating slabs of dimen- 
sions 130 X 2.5 cm^, with a thickness of 2.5 cm. Both ends of each slab are 
coupled to specially shaped light-guides as illustrated in Fig. 24. 

Photoconversion is done by Philips XP2008 photomultipliers, equipped, as for 
the VH, with a voltage divider allowing high counting rate capability. The 
front-end electronics system contains the same elements as those described for 
the vertical hodoscopes. The single hit detection efficiency of HH is greater 
than 96.6 % on both arms, and the time resolution is 320 ps. 



11 The Cherenkov Counters 



This detector is essential for rejection of the main background of electron- 
positron pairs from photon conversion, Dalitz pairs, and to a minor extent 
from resonance decays. It is used in the first level tt+tt" main trigger and in 
the calibration trigger to select e'^e~ pairs. 
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It is structured in two identical threshold Cherenkov counters [28], each cov- 
ering one spectrometer arm (see Fig. 25). 



The gas radiator is enclosed in a volume defined by the entrance and exit 
windows, with dimensions 143 x 56 cm^ and 336 x 96 cm^, respectively. The 
chosen radiator is N2 at normal temperature and pressure {9(^ — 1.4°) and the 
counter length is 285 cm. 




Fig. 25. The far end part of the DIRAC setup, comprising threshold Cherenkov 
counters (Ch), preshower detector (PSh), iron absorber (Fe) and muon counters 
(Mu). 

Each counter is equipped with 20 mirrors and 10 photomultipliers on two 
rows. Cherenkov light reflected by pairs of adjacent mirrors is focused onto the 
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Fig. 26. ADC spectrum from one Cherenkov photomultiplier: (a) spectrum from 
pions (practicahy equal to the ADC pedestal distribution), (b) amphtude signal 
from single photoelectron, (c) spectrum from electrons. 
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Fig. 27. Distribution of the number of photoelectrons detected from the (a) positive 
and (b) negative Cherenkov detector arms. 

same photomultiplier (Hamamatsu R1587, with 130 mm UV-glass window). 
Mirrors and PM arrays are placed along circumferences with radii 630 cm 
and 570 cm, respectively, centred at the magnet middle point. Mirrors are 
spherically deformed rectangles, with average dimensions 30 x 35 cm^ and 
6 mm thickness. 

The analog signals from individual PM are fed into two custom-made summing 
modules, one per counter (10 channels input). The output of the summing 
module is a linear sum with fan-out of 2, plus the 10 input signals. At the 
module output the linear analog sum is attenuated by a factor of 3 with respect 
to individual channels. One analog sum is discriminated (LeCroy 4413 leading 
edge discriminators) and used for trigger purposes, whereas the individual 
PM signals and the other sum are fed, after a delay line, into LeCroy 4300B 
ADC units. Individual PM ADC spectra are used to perform overall amphtude 
alignment on the basis of maximum ratio between mean and r.m.s., in the PM 
voltage range 2400 ± 200 V. 

A single photoelectron peak is clearly observed in all channels and it is shown 
in Fig. 26(b). The single photoelectron spectrum was used to cross-check 
the conversion factor from ADC counts to number of photoelectrons {N^g), 
which is obtained from the analysis of the summed amplitude spectra. Fig. 27 
shows the latter as a function of Np^ for the positive (a) and the negative 
(b) spectrometer arms. The mean values are A'pe=16.2 and A'pe=16.4, respec- 
tively. From these values we infer that both counters have an efficiency greater 
than 99.8 % when operated at a threshold shghtly less than 2 photoelectrons. 
The pion contamination above the detection threshold is estimated to be less 
than 1.5 %. Such contamination arises from pions with momenta above the 
Cherenkov threshold and from accidental coincidences occurring within the 
trigger time- window. 
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Fig. 28. General picture of the spectrometer installed downstream the magnet, taken 
from a reference plane close to the Cherenkov counters, which can be appreciated 
in the bottom part of the picture. The magnet can be seen at the far end, and in 
between we see the vertical and horizontal hodoscope photomultiplier layout. The 
drift chamber system is installed behind it. 

A general picture of the installed detectors described so far, looking into the 
magnet from a distance close to that of the Cherenkov counter, can be appre- 
ciated in Figure 28. 



12 The Preshower Detector 

The purpose of the Preshower detector (PSH) is two-fold: it provides addi- 
tional electron/pion separation power in the off-line analysis and is used in 
the trigger generating logic as well. 




Fig. 29. PSH element and PSH array on one spectrometer arm. 

The PSH is based on an array of lead converters followed by scintillation 
detectors [29]. Electrons (positrons) initiate in the converters electromagnetic 
showers which are sampled in the scintillation counters while pions behave 
mainly as minimum ionising particles. For trigger purposes, the signal has to 
be produced whether a pion or an electron cross a PSH counter. 

The PSH consists of 16 detector elements placed symmetrically in two arms, 
as seen in Fig. 25. Each element has a Pb converter and a scintillation counter. 
Fig. 29. The converters of the two outermost elements of each arm (low mo- 
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mcntum region) arc 10 mm thick, whereas the rest are 25 mm thick (around 
2 and 5 units of radiation length, respectively). The scintillator used is BI- 
CRON type BC-408, with slab dimensions 35 x 75 cm^ and 1 cm thickness. 
The scintillation hght is transmitted to photomultiphers EMI 9954-B, placed 
at one end only, by 10 mm thick Plexiglas light-guides ending with Plexiglas 
cylinders to match with the PM photocathodes. Since the maximum particle 
flux on each PSH clement is as high as 2 MHz, an additional booster power 
supply is used to feed the last PM dynodes. 

The detector signals are hnearly split into two branches, one used for trigger 
purposes and another for ADC analysis. In the former, a leading edge LcCroy 
4416 discriminator is used with a threshold corresponding to efficient detection 
of minimum ionising particles. In the latter the signals are fed into 20 dB 
attenuators in order to fit the electron signals into the dynamical range of the 
LeCroy 4300B ADC. 

The single arm detector efficiency is 99.5% for pions. Fig. 30 shows the pulse- 
height spectra from one element of the PSH for pions and electrons selected by 
the trigger system (this selection is based on the Cherenkov detector response) . 
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Fig. 30. Pulse-height spectra for pions and electrons in one element of PSH. 



As seen from Fig. 30, the pion spectrum has a tail (originating from nuclear 
interaction of pions in the converter) extending to the electron amplitude 
region. The off-line study of the c/tt separation efficiency of the PSH showed 
that rejection of electrons reaches 85% with less than 5% losses of pions. As 
only a very small fraction of electrons escapes the on-line identiflcation in the 
Cherenkov counters, the combined use of the Cherenkov and PSH detectors 
provides almost 100% electron rejection power at the off-hne analysis stage. 
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13 Muon detector 



Admixture of muons in the tttt events can be a serious source of background. 
For this reason a muon detection system is implemented to provide efficient 
muon tagging. Muons come almost entirely from pion decays with a small 
admixture from other decays and direct ii'^/J," pair production. 

The muon detector consists of scintillation counters placed behind a thick iron 
absorber which almost entirely absorbs hadrons and related hadronic showers. 
This detector is placed at the downstream end of the DIRAC apparatus, few 
meters from the intense primary proton beam dump. As a result, the muon 
scintillation counters may undergo a high flux of background radiation from 
the beam dump area. This has required a special design of the counter arrays 
and electronics and has prevented from using muon information during on-line 
data reduction. 



The counters are located behind iron absorber blocks with thickness ranging 
from 60 to 140 cm (see Fig. 25). The thickness is larger in the region close 
to the spectrometer symmetry axis, in order to compensate for the harder 
pion momentum spectrum. A double layer structure has been envisaged for 
the counters, each layer on each arm consisting of 28 counters with equal 
scintillating slabs of 75 x 12 cm^ front area and 0.5 cm thickness. The muon 
detector data are read out only if simultaneous signals from a pair of corre- 
sponding counters in the two layers are detected. This essentially reduces the 
background counting rate induced by the neutron flux from the beam dump. 
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Fig. 31. Readout scheme of the muon detector element. 
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Fig. 32. Schematic layout of muon counters on their support structure, indicating 
hght guides and photomultiphers. 

A special readout architecture is realised to compromise between a cost-saving 
solution and the need to achieve a reasonably high time resolution. Scintilla- 
tion light is detected by 25 mm diameter bialkaline FEU-85 photomultiphers 
placed at one scintillator end, in the two layers at opposite ends, as shown in 
Fig. 31. Signals from a pair of counters are fed into constant fraction discrim- 
inators (CFD) followed by meantimer (MT). CAEN modules C808 and C561 
are used, respectively, for this purpose. In such a scheme the output signal is 
generated only if both counters arc hit, and correct timing occurs only if the 
same particle crosses both counters. 

Fish-tail light guides are used to couple the PM photocathodes to the scintil- 
lators, except when this is impeded by the presence of the concrete floor (see 
Fig. 32). 

In this case (for 20 counters of the second layer in each arm), photomultiph- 
ers are directly coupled to the scintillating slabs, which are then made twice 
thicker to compensate for the loss of light yield. 

In Fig. 33 the measured time difference between the signals of the muon de- 
tector and the vertical hodoscope is shown for an individual muon element of 
the positive arm. A global time resolution of 1.3 ns is observed, with a very 
small background level. 

The way to suppress muons at the trigger level would be to include the muon 
counter signals into the anti-coincidence trigger logic. Because of the high 
background load to the muon detector, one might dangerously suppress useful 
pion events if they happen to occur on-time with background signals in the 
muon detector. That is why the option to use off-line the muon detector in- 
formation has been chosen. In the off-line analysis only the events with muon 
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Fig. 33. Time difference between tfie signals of the muon detector and the vertical 
hodoscope. 

counter hits correlated in time with those of other detectors are tagged as 
"muon" events and thus rejected [30]. From the analysis of experimental data 
we have inferred that the fraction of such events, containing at least one muon, 
is about 10% [31]. Such muon-events originate to a large extent (~ 80%) from 
TT^-decays in the path between the DC and the muon counter systems (de- 
cays upstream the DC are mostly suppressed by the trigger system and thus 
contribute less to the collected event sample). 



14 Trigger system 

The trigger system was designed to provide a reduction of the event rate 
to a level acceptable to the data acquisition system which is around 2000 
events/spill. Pion pairs are produced in the target mainly in a free state with 
a wide distribution over their relative momentum Q, whereas atomic pairs 
from disintegration have very low Q, typically below 3 MeV/c. The on- 
line data selection rejects events with pion pairs having > 30 MeV/c or 
Qx > 3 MeV/c or Qy > 10 MeV/c, keeping at the same time high efficiency 
for detection of pairs with Q components below these values, {Ql,Qx and 
Qy are longitudinal and transversal components of the relative momentum, 
respectively.) 
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A multilevel trigger is used in DIRAC [32] . It comprises a simple and fast first 
level trigger and higher level trigger processors which apply selection criteria 
to different components of the relative momentum of pion pairs. 

Due to the requirements of the data analysis procedure, the on-line selection 
of only time correlated (prompt) pion pairs, originating from a single proton- 
target interaction and detected simultaneously by both spectrometer arms, 
is not enough. In addition, a large number of uncorrelated, accidental, pion 
pairs is also necessary. These accidental pairs are used in the off-line analy- 
sis to describe the relative momentum distribution of free (non-atomic) pion 
pairs without Coulomb interaction in the final state. Therefore, the trigger 
system should apply very similar selection criteria to prompt and acciden- 
tal events, within a preselected coincidence time window centred around the 
peak of prompt events. The statistical error of the A2Tr lifetime measurement 
depends on the number of both prompt and accidental detected pairs. In 
standard experimental conditions, the optimal ratio of prompt to accidental 
events is obtained using a 40 ns wide coincidence time window between the 
times measured in the left (VHl) and right (VH2) vertical hodoscopes. 

Since 1999, when the experiment has started, the trigger architecture was 
upgraded several times to achieve a larger reduction of the background event 
rate (prompt and accidental pairs with large values of Q). In the present article 
we briefly describe the most recent version ^ . 

A block diagram of the trigger architecture is presented in Fig. 34. The first 
level trigger Tl starts digitisation of the detector signals in the data acquisi- 
tion (DAQ) modules (ADC, TDC, etc.). At the next level the neural network 
trigger DNA/RNA (DIRAC Neural Atomic and Revised Neural Atomic trig- 
ger) rejects the events with high Q values. At the last stage, a powerful drift 
chamber trigger processor T4 imposes additional constraints to the relative 
momentum and takes the final decision to accept or to reject the event. 

In addition to the main trigger designed to detect pionic atoms, several cali- 
bration triggers are run in parallel. They are applied to DAQ directly, without 
suppression by higher level trigger stages but with appropriate prescaling fac- 
tors. 



Before 2001 the trigger system included T2 and T3 stages following the first level 
trigger (Tl). T2 selected particle pairs with a small Ax distance in the upstream 
spectrometer region (for rejection of high Q^:) using the data from SFD and IH. T3 
analysed hit patterns of the upstream IH and downstream VH detectors imposing 
selection criteria to Ql values. A detailed description of T2 and T3 is given in [32] 
and [33] . With the implementation of a new architecture of the trigger system these 
trigger stages were removed. 
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Fig. 34. General block diagram of the DIRAC multilevel trigger. 

14.1 First level trigger (Tl) 



The first level trigger (a detailed description is given in [34]) fulfils tfie following 
tasks: 

— Selects events with signals in both detector arms downstream the magnet. 

— Classifies the particle in each arm as or depending on the presence 
of the Cherenkov counter signal. Protons, kaons and muons are equally 
included in the "pion" class, their identification is performed in the off-line 
analysis. 

— Arranges the coincidences between the signals detected in the two arms. 
The width of the coincidence time window defines the ratio between yields 
of prompt and accidental events in the collected data. 

— Applies a coplanarity criterion to particle pairs: the difference between the 
hit slab numbers in the horizontal hodoscopes in the two arms (HHl and 
HH2) should be < 2. This criterion forces a selection on the Qy component 
of the relative momentum and provides a rate reduction by a factor of 2. 

— Selects in parallel events from several physics processes needed for the setup 
cahbration: e+e" pairs, A — > p + tt" decays, decays to three charged 
pions. 



The physics and calibration trigger signals pass through the mask register 
and, after proper prescaling, are combined with an OR function. Any trigger 
type can be enabled or disabled by proper programming of the mask register. 
Independent prescaling of each sub-trigger channel allows to adjust their rel- 
ative rate with respect to the rate of the main trigger. A specific trigger mark 
is recorded for every event to allow sorting the data by trigger type during 
off-line analysis and on-line monitoring. 

All Tl modules are ECL line programmable multichannel CAMAC units. Most 
of them are commercial modules, except for the dedicated coplanarity proces- 
sor which has been custom-developed at JINR. Meantimer units are used in 
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all VH and HH channels to remove the dependence of the time measurement 
on the hit location, thus reducing the total trigger time jitter. 



14.2 Neural network trigger (DNA/RNA) 



The DNA/RNA trigger [35] is a processing system using a neural network 
algorithm. Its hardware is based on the custom-built version of the neural 
trigger used in the CPLEAR experiment [36]. 

DNA/RNA receives (see Fig. 35) the hit patterns from the vertical hodoscopes 
VHl, VH2 and the X-planes of the upstream detectors: the ionisation ho- 
doscope (IH) and the scintillating fibre detector (SFD). For low Q events the 
hits in these detectors are correlated. 




to RNA trigger 

Fig. 35. DIRAC detectors used for the neural network triggers DNA and RNA. 
Numbers of signal channels from each detector are given in parentheses. 
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The neural network was trained to select particle pairs with low relative mo- 
menta: < 3 MeV/c, Qj^ < 10 MeV/c and Ql < 30 McV/c. The events 
which do not satisfy any of those conditions are considered "bad" and re- 
jected. 

The DNA/RNA logic is started by a fast pretrigger, TO, and in 250 ns evaluates 
an event. Rate reduction by a factor of 2 with respect to Tl is achieved with 
DNA/RNA. 



14.3 Drift chamber processor (T4) 



T4 is the final trigger stage. T4 processor reconstructs straight tracks in the 
X-projection of the drift chambers and analyses them to determine the value 
of the relative momentum (the algorithm is described in [32] ) . 




Fig. 36. T4 operation block diagram. Only the drift chamber X-planes involved in 
T4 are shown. 



The block diagram of the T4 operation is shown in Fig. 36. The drift chamber 
processor includes two stages: the track finder and the track analyser. The 
track finder (an identical processor is used for each arm) receives the numbers 
of the hit wires from all drift chamber X-planes. Drift time values are not 
used in the T4 logic. A unique number, "track identifier" , which contains the 
encoded numbers of the hit wires, is associated to the found track. 
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If tracks arc found in both arms, the track analyser continues the event eval- 
uation. The track analyser receives the track identifiers from both arms and 
compares them with the content of a look-up memory table which contains all 
possible combinations of track identifiers for pion pairs with Ql < 30 MeV/c 
and Qx < 3 MeV/c. These "allowed" combinations are obtained from a ded- 
icated simulation using the precise geometry of the setup. If a relevant com- 
bination is found, the T4 processor generates a positive decision signal which 
starts the data transfer to the VME buffer memories. Otherwise, the Clear 
and Reset signals are apphed to the DAQ and trigger systems. 

The T4 decision time depends on the complexity of the event and is around 
3.5 /is on average. The rejection factor of T4 is around 5 with respect to the 
Tl rate and around 2.5 with respect to DNA/RNA. 

The whole trigger system is fully computer controlled: no hardware interven- 
tion is needed in order to modify the trigger configuration. With all selection 
stages enabled the event rate at the typical experimental conditions is around 
700 per spill, that is well below the limits of the DAQ rate capability. 



15 Data acquisition system 



The architecture of hardware and software parts of the data acquisition system 
takes into account the time structure of the proton beam. The machine super- 
cycle of the CERN PS accelerator has 15-20 s duration. Within this period 
DIRAC receives protons in spills of 400-450 ms width, from 1 to 5 spills per 
super-cycle. The interval between consecutive spills can be as short as 1 s. 

During the accelerator burst the data from all detectors are read out into 
VME buffer memories (commercial and dedicated electronic modules) with- 
out any software intervention. The data transfer to VME processor boards, 
event building, data transfer to the main host computer and other relatively 
slow operations are performed during the pause between bursts. This provides 
the maximum operation rate of DAQ. The information comes from 2048 chan- 
nels of microstrip gas chambers, 800 channels of scintillating fibre detector, 
2016 channels of drift chambers and 224 channels of other scintillation and 
Cherenkov detectors. For every channel the time or amplitude information is 
recorded or both. In addition to a main readout mode, the readout of scalers 
at the end of every spill is arranged via a CAMAC bus. The counting rates 
of all the detectors, trigger rates at different trigger levels and for different 
sub-trigger modes are recorded together with the value of the beam intensity 
provided by the PS complex. 
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15.1 DAQ hardware 



The data readout [37] is arranged with 12 readout branches: 4 branches for 
GEM/MSGC, 3 branches for DC and 5 FERA [38] branches for all other 
detectors. In FERA and DC branches the VME modules CES HSM 1870 and 
LeCroy 1190 are used as buffer memories. In GEM/MSGC branches the buffer 
memories are incorporated into dedicated VME modules [39] . FERA branches 
include different FERA compatible LeCroy modules: ADC4300B, multi-hit 
TDC3377, universal logic modules 2366 configured in this application like 
FERA registers and scalers. Pecuharities of FERA readout in DIRAC, such 
as multi-gate and Fast Clear operation, are described in [40] . 

The logic of readout is the following. The first level trigger Tl starts digiti- 
sation in ADC and TDC of FERA and DC branches (see Fig. 34 in Section 
"Trigger system" ) . Readout is inhibited unless a positive decision of the high- 
est level trigger T4 is received. If the T4 processor decides positively, the 
Inhibit Readout status is released and the converted event data arc trans- 
ferred to buffer memories. If the decision of DNA/RNA or T4 is negative, 
then a Fast Clear signal is generated [40] which discards the data in all FERA 
modules and DC registers. 

In contrast to FERA and DC subsystems, the processing of GEM/MSGC data 
is started by the next level DNA/RNA trigger. This reduces the dead time in- 
troduced by the Clear process which in GEM/MSGC electronics takes longer. 
Thus, the negative decisions of T4 only lead to clearing of the GEM/MSGC 
data. If no Clear signal is received, the converted data are transferred to buffer 
memories. 

The readout of the whole event takes 45 /is and is defined by a fixed acquisition 
time of GEM/MSGC which exceeds the readout time in other branches. 



15.2 DAQ software 

The main part of the DAQ software [41] is running on two VME processor 
boards and on the main DAQ host. PowerPC-based VME processor boards 
control the trigger and FERA readout electronics via two CAMAC branch 
drivers, operate VME modules, read data from buffer memories and transfer 
them to the main DAQ host. The main DAQ host performs event building, 
records the built data and distributes them to other computers for on-line 
monitoring and analysis. 
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The DAQ software is written in C programming language and is running under 
UNIX-like operating systems: Lynx-OS (VME processors) and LINUX (main 
DAQ host). 

The schematic layout of the data acquisition processes is presented in Fig. 37. 
The basic processes running on VME and main DAQ hosts and on one of 
the monitoring computers are shown. Only one VME processor is shown for 
simplicity. 



DAQ host 




Fig. 37. Layout of the DAQ architecture. 

The DAQ software can be divided into more or less independent groups. The 
first one is a set of programs for data readout (runs on VME boards) and 
event building (runs on the main DAQ host). During the event building the 
data are checked for consistency using the headers and the serial numbers of 
sub-events provided by the detector readout electronics. The serial numbers 
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allow to check that all buffer memories contain the same number of sub-events 
and that different sub-events correspond to each other. 

The second group consists of programs for data distribution over the net- 
work for on-line monitoring. The group includes a primary data receiver and 

secondary data senders/receivers. The primary receiver is the only program 
which receives data directly from the event builder and is critical for the data 
acquisition: other processes of this group may either be or not be running and 
may be started or stopped at any time. 

The third group of programs allows to control the run status via a graphical 
user interface: to set the run parameters, to select the trigger type, to change 
the run status (start, stop, suspend etc.). This group includes also a Run 
Display program which delivers information about the current run conditions. 

The fourth group is a set of slow-control processes providing a uniform access 
to hardware from any host participating in DAQ. The tasks of these programs 
are, for example, high voltage control, VME and CAMAC module initialisation 
and so on. 

The whole DAQ system is able to accept data from five consecutive bursts in 
one accelerator super-cycle with up to 2 MBytes per burst (this corresponds to 
2000 triggers per burst) and to distribute events to all hosts participating in 
on-hne data processing. The limitation comes from the capacity of the VME 
buffer memories. 



15.3 On-line monitoring 

An on-line monitoring program receives the data distributed by the main 
DAQ process. The program fulfils accumulation of several hundred histograms, 
including all raw signal spectra from any detector as well as spectra obtained 
after a fast preliminary data handling. The software is written in the frame of 
ROOT software package. 

The program delivers the following information: 

— Time and amplitude distributions for every counter. 

— Hit and multiplicity distributions. 

— Correlation plots for hits in different detectors. 

— Beam profile on the experimental target, beam spill time structure. 

— Some specific data for selected detectors (cluster size for GEM/MSGC etc.). 

— Event rates for different trigger modes running in parallel and their timing. 
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— Information about the data taking process as a whole: burst data volume, 
event volume, time of event in burst etc. 

Proper setting of the corresponding program parameter allows selection of the 
trigger type for which all the above operations should be fulfilled. For example, 
the histograms can be obtained separately for e'^e~ , tt'^tt", A and any other 
trigger or for all accepted triggers. 

Apart from the common program, there are dedicated monitoring programs 
developed for more detailed control of individual detectors. 



16 Setup performances 

The experimental setup is designed to provide high efficiency for detection of 
7r"'"7r~ pairs with small opening angle in the laboratory frame {6 < 3 mrad) 
and small relative momentum in the center of mass frame {Q < 3 MeV/c). 
Detector resolution and multiple Coulomb scattering in the target and setup 
elements affect the measurement accuracy. 

Detector resolution leads to a linear increase of the relative error (Tp/p with in- 
creasing momentum value, whereas multiple scattering contributes a constant 
term to the momentum uncertainty. With an experimental space resolution of 
90 ^m in the DC system and 50 ^m in the GEM/MSGC detectors, the uncer- 
tainty on the coordinate and momentum measurements are mainly determined 
by the description of the multiple scattering errors. 

Both momentum and the opening angle measurement accuracies affect the 
precision on the measurement of the relative momentum Q of the pion pair. 
The momentum accuracy ap/p is 0.3%, at Pt^ — PA2n/'^ — 2 GeV/c, almost 
independent on momentum. 

The kinematics of tt'^tt" pairs originating from the ionisation of A2Tr atoms in 
the target requires ability on two-track separation in the upstream detector 
region better than 1.5 mm in most of the cases. The experimental double track 
resolution is limited by the granularity and signal clustering of the upstream 
detectors and with current algorithms is about 0.4 mm. Thus, the accuracy 
on the measurement of the longitudinal (Ql) and transversal components {Q^ 
and Qy) of Q is 0.6 MeV/c (Ql) and 0.4 MeV/c {Q^ and Qy). These errors 
arise from the estimated contribution of multiple scattering in the detectors 
and setup elements, according to our present knowledge of the experimental 
data. In fact, the effect of multiple scattering in the foil target increases the 
error in and Qy to up 1 MeV/c. 
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Fig. 38. Invariant mass distribution of pvr pairs. Events were selected from a sample 
of A calibration data. 



The calibration of the energy scale and the experimental determination of 
the setup resolution are performed by monitoring the position and the width 
of the A detected by the apparatus. The reconstructed p7i~ invariant mass 
obtained from calibration data is shown in Fig. 38. The A is clearly seen 
above a small background. A gaussian plus a polynomial fit to the data gives 
Ma = 1115.67 MeV/c^ and cta = 0.43 MeV/c^. Calibration using A pir' 
decays is used in addition to control the precision of the setup alignment. 
Any misalignment of the tracking system in one arm relative to the other arm 
would result in asymmetrical errors on the reconstructed momenta. 
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Fig. 39. Left: Distribution of Q for accepted vr^vr^ pairs after the full DIRAC trigger 
system, and for minimum bias pairs. Right: trigger acceptance, determined as the 
ratio between the previous two distributions. 
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Fig. 40. Experimental p and pT distributions of vr corrected for apparatus accep- 
tance (histogram) with superimposed the results of the parameterisation [42] [43] . 

The accuracy of the time measurement is obtained from the analysis of the 
distribution of the e^e~ time difference at the vertical hodoscopes as explained 
in that section (see Fig. 22). 

The performance of the trigger system as a whole in selecting low-Q events 
is illustrated in Fig 39, where the magnitude of relative momentum of pion 
pairs Q is shown (in their center-of-mass frame), after DNA/RNA and T4 
trigger selection (see trigger section). Trigger efficiency as function of Q is flat 
in the low-Q region, as illustrated in Fig. 39. This is considered an important 
figure of merit of the spectrometer, for a precision study of the tt"*"?!" Coulomb 
interaction. 

Pions from ionisation of A2t^ entering the apparatus have momenta below 
4 GeV/c. The apparatus momentum acceptance for time-correlated pairs is 
fiat for pions with momenta between 1.6 GeV/c and 3 GeV/c, and it decreases 
for higher momenta. 

For the sake of completeness we also show the p and Pt distributions for a 
single TT^ in Fig. 40, with superimposed a parameterisation of the inclusive 
yield based on the analytic representation [42], adapted to DIRAC center-of- 
mass energy (^=6.84 GeV) [43]. 
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